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In the past decade imaging atmospheric Cherenkov telescope arrays such as H.E.S.S., MAGIC, 
VERITAS, as well as the Fermi-LKY space telescope have provided high quality images and 
spectra of the 7 -ray universe. Currently the 7 -ray community is preparing to build the next- 
generation Cherenkov Telecope Array (CTA), which will be operated as an open observatory. 
Gammapy vO.3 (available at https://github.com/gammapy/gammapy under the open-source BSD 
license) is a new in-development Astropy affiliated package for high-level analysis and simulation 
of astronomical 7 -ray data. It is built on the scientific Python stack (Numpy, Scipy, matplotlib and 
scikit-image) and makes use of other open-source astronomy packages such as Astropy, Sherpa 
and Naima to provide a fiexible set of tools for 7 -ray astronomers. 

We present an overview of the Gammapy scope, development workfiow, status, structure, features, 
application examples and goals. We would like Gammapy to become a community-developed 
project and a place of collaboration between scientists interested in 7 -ray astronomy with Python. 
Contributions welcome! 
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Figure 1: Gammapy is a Python package for high-level y-ray data analysis. Using event lists, exposures 
and point spread functions as input you can use it to generate science results such as images, spectra, light 
curves or source catalogs. So far it has been used to simulate and analyse H.E.S.S., CTA and Fermi-LP^ 
data, hopefully it will also be applied to e.g. VERITAS, MAGIC or HAWC data in the future. 

1. Introduction 

1.1 What is Gammapy ? 

Gammapy is an open-source Python package for y-ray astronomy. Originally Gammapy started 
as a place to share morphology fitting Python scripts for the work on the H.E.S.S. Galactic plane 
survey [1] two years ago. Since that time Gammapy has grown steadily: functionality as well 
as development infrastructure has been improved and it has been accepted as an in-development 
Astropy-affiliated package. Now, in this proceeding, we would like to introduce Gammapy to the 
community and present our vision of Gammapy as a future community-developed, general purpose 
analysis toolbox for y-ray astronomers. 

The general concept of Gammapy is illustrated in Figure 1. Based on pre-processed input 
data (e.g. event lists) provided by instruments such as H.E.S.S., Fermi or CTA, Gammapy offers 
the high-level analysis tools to generate science results such as images, spectra, light curves and 
source catalogs. By using common data structures and restriction to binned analysis techniques, all 
input data can be treated the same way, independent of the instrument. Research already making 
use of Gammapy is presented in [1, 2, 3]. 

1.2 How to get Gammapy 

Recently Gammapy version 0.3 was released. It is available via the Python package index ^ or 
using package manager tools like pip and conda. Gammapy works on Linux and Mac (Windows 
most likely as well, but this has not been tested yet) and is compatible with Python 2.7 and Python 

^ https ://pypi.py thon.org/pypi/gammapy/ 
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Figure 2: The Gammapy stack. Required dependencies Numpy and Astropy are illustrated with solid 
arrows, optional dependencies (the rest) with dashed arrows. 


3.3 or later. Further details on requirements and installation are available online The latest 
documentation is available on Read the Docs including tutorials, code and analysis examples. 
We have also set up a mailing list for user support and discussion 


2. The Gammapy stack 

Gammapy is primarily built on the scientific Python stack. We employ Numpy and Astropy [4] 
as main dependencies and integrate other packages as optional dependencies where necessary. The 
current dependency structure is illustrated in Figure 2. Numpy provides the low level data struc¬ 
tures and the framework for numerical/array computations. Astropy is used for higher level data 
structures like tables (Table) and n-dimensional data objects (NDData), for I/O, coordinates and 
WCS transformations, and handling of physical quantities with units. Scipy is used for advanced 
numerical and data processing algorithms. For specific image processing routines we additionally 
use scikit-image. 

To allow morphology and spectral fitting of TeV sources with Gammapy we use the established 
X-ray modelling and fitting package Sherpa [5]. Sherpa allows interactive and scripted fitting of 
data sets with various spectral, light curve and morphology models taking instrument response 
functions into account. Additionally, it is possible to determine confidence levels on best-fit model 
parameters, compute likelihood profiles and goodness of fit measures. Sherpa recently also became 
an open-source project which makes it possible for users and external developers to contribute 
missing functionality or fix issues themselves in future. 

^ http s ://gammapy.readthedoc s. org/en/latest/install.html 
^ http s ://gammapy.readthedocs. org/en/latest 
https ://groups. google, com/forum/#! forum/gammapy 
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Modelling and fitting of non-thermal radiation processes to spectral energy distributions (SEDs) 
is provided with the Naima package. It uses Markov-Chain Monte Carlo emcee sampling to find 
best-fit model parameters of physical radiation models and thus determine the radiation mecha¬ 
nism that leads to the observed emission. The radiative models available in Naima can be called 
as source models in the Gammapy environment and be used, for example, to simulate spectra for 
source population studies. 

Data visualizing and plotting of sky images is done with matplotlib and the astropy-affiliated 
package wcsaxes. Further astropy-affilated packages we allow as optional dependencies are pho- 
tutils and reproject. Photutils is used for source detection and photometry and reproject for re¬ 
projection of sky images. 

While this is a large number of dependencies, most of them are optional and only needed for 
small specific tasks so we do not see this as a strong disadvantage. Optional packages can be easily 
installed using package managing tools like pip or conda. By re-using other packages’ functionality 
which would be complex to re-implement, new tools and techniques can be easily integrated into 
Gammapy to help obtain scientific results more quickly with a minimal coding effort. 

3. Development workflow 

Gammapy uses all the standard tools of modern community-driven open-source software de¬ 
velopment. The code repository is hosted on GitHub which allows for convenient and public 
interaction of developers, contributors and users via the GitHub interface. This includes, for in¬ 
stance, requests for missing functionality, issues for bugs or questions and pull requests for 
code contributions. The code base is continuously built and tested in different virtual environments 
using the continuous integration service Travis CL This helps to maintain high coding standards 
and compatibility with different versions of Numpy and Astropy. 

The documentation of Gammapy is generated using Sphinx. An online version of the latest 
documentation is built automatically and is available on Read the Docs In addition to the main 
code repository we maintain gammapy-extra. This repository contains prepared catalog and map 
data and IPython notebooks with more extensive analysis examples and tutorials. 

4. The Gammapy toolbox 
4.1 Sub-packages 

The Gammapy code base is structured into several sub-packages where each of the packages 
bundle corresponding functionality in a namespace. This follows the concept of other Python 
packages such as Astropy and Scipy. The following list gives a rough overview of the different 
sub-packages with a short description: 

• gammapy. astro - Galactic population and emission models of TeV sources 

• gammapy. background - Background estimation and modeling 

^ http s ://github. com/gammapy/gammapy 

^ http s ://gammapy.readthedocs. org/en/latest 


4 




Gammapy - A Python package for y-ray astronomy 


Axel Donath 


• gammapy. catalog - y-ray source catalog access and processing 

• gammapy. datasets - Easy access to bundled and remote datasets 

• gammapy. detect - Source detection tools and algorithms 

• gammapy. hspec - Interface to spectral fitting with Sherpa 

• gammapy. image - Image processing and analysis tools 

• gammapy. irf - Instrument response function (IRE) access and handling 

• gammapy .morphology - Morphology models and tools 

• gammapy. obs - Observation handling 

• gammapy. spectrum - Spectrum models and tools 

• gammapy. stats - Statistical functions 

• gammapy. time - Handling of time series and y-ray lightcurves 

• gammapy. utils - Utility functions and classes (in sub-modules) 

There are some cases (e.g. gammapy. spectrum, models and several sub-modules of 
gammapy. utils) where the end-user functionality is exposed one level further down in the hierar¬ 
chy. This is because putting everything into the top-level gammapy. spectrum or gammapy .utils 
namespace would lead to an unstructured collection of functions and classes. A large part of 
Gammapy ’s functionality uses an object-oriented API. Eunctions are only used where it leads 
to an improved or more intuitve user interface. 

4.2 Command-line tools 

Gammapy includes various ready to use command-line tools which provide a familiar interface 
to data processing for many astronomers. A complete overview of all available tools can be found 
online ^. When Gammapy is installed the user can simply type gammapy- on the command-line 
and use tab completion to see the list of tools, as all Gammapy tools start with the same gammapy- 
prefix. Specific information on single command-line tools and a description of available parameters 
is shown when calling the corresponding tool with the standard -h or — help option. 

5. Usage examples 

The Gammapy package should be considered as a toolbox out of which powerful analysis 
scripts can be composed easily by astronomers, even if they have very little programming experi¬ 
ence. In the following section we present a selection of code examples demonstrating how to set 
up complex analysis steps with just a few lines of code. 

5.1 Morphology fitting command-line tool 

Gammapy includes a simple to use command-line script gammapy-sherpa-like to perform 
a Poisson maximum likelihood morphology fitting using EITS files as input. The input counts, 
exposure and background maps can be specified by the corresponding parameters — counts, 
— exposure and — background. The source model is defined in a ISON file and should be 
passed using the — sources option. The model parameters for a multi-Gaussian point spread 
function (PSE) can be specified in ISON format using the — psf option. 

^https://gammapy.readthedocs.org/en/v0.3/scripts/index.html 
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Figure 3: Fermi survey TS map. 


1 $ gammapy-sherpa-like —counts counts.fits —exposure exposure.fits 

2 —background background.fits —psf psf.json —sources sources.json 

3 result.json 


Fit results and additional information are stored in result. json. 


5.2 Test statistics maps 

The gammapy. detect module includes a high performance compute_ts_map function to 
compute test statistic (TS) maps for y-ray survey data. The implementation is based on the method 
described in [6]. As input data, the user provides counts, background and exposure maps. The 
following code example demonstrates the computation of a TS map for prepared Fermi survey 
data, which is provided in ganunapy-extra. The resulting TS map is shown in Figure 3. 


1 from astropy.io import fits 

2 from astropy.convolution import Gaussian2DKernel 

3 from gammapy.detect import compute_ts_map 

4 hdu_list = fits.open( 'all.fits.gz' ) 

5 kernel = Gaussian2DKernel( 2.5 ) 

6 result = compute_ts_map(hdu_list[ 'On' ]•data^ 

7 hdu_list [' Background' ].data, 

8 hdu_list[ 'ExpGammaMap' ].data, kernel) 

9 result.write (' ts_map.fits' ) 


5.3 Galactic population models 

Gammapy also offers functionality for modeling and simulation. The gammapy .astro sub¬ 
package provides tools to simulate Galactic TeV source populations and source characteristics, 
which is useful in the context of surveys and population studies. Figure 4 illustrates a source pop¬ 
ulation simulated with Gammapy and the different radial distribution models the user can choose 
from. 
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Figure 4: Galactic source population simulated using 
after [7] and the spiral-arm model of [8]. 
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Gammapy , assuming a radial distribution of sources 


6. Planned functionality 

Gammapy already includes key analysis features like morphology and spectrum fitting but 
both are currently limited to either 2D image-based data or ID spectral data. As a major next step 
we plan to support joint likelihood fitting of datasets, as illustrated in Figure 5. Events are binned 
into longitude, latitude and energy cubes and fitted simultaneously with spectral and spatial models 
taking energy-dependent background, exposure and point spread function (PSF) into account. 

This approach will allow joint likelihood analysis across different experiments. For instance 
simultaneous likelihood fitting of Fermi and H.E.S.S. data. 

Support for un-binned analyses is not planned. 


7. Summary 

Gammapy 0.3 is still alpha quality software. It is a package where standard y-ray analyses are 
available on the one hand while, on the other hand, integration and prototyping of new methods is 
easily possible. So far it only contains limited functionality but the setup of documentation, testing 
and deployment is already very advanced. It’s scope will continously grow and we hope that 
many users and developers show interest in open and reproducible y-ray astronomy with Python. 
As long-term goal we would like Gammapy to turn into a fully community-developed package. 
So all contributions to Gammapy are welcome! If you don’t know how to turn your scripts into 
production-quality, reusable code, please get in touch with us (e.g. using the mailing list) and we 
will help you get there! 

For further information on how Gammapy and other tools are being used for H.E.S.S. data 
analysis, we encourage you to look at [9]. 
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Figure 5: Gammapy data model illustration. Binned analysis of lon-lat-energy cube data is supported via 

joint likelihood analysis of one image per energy bin. On-off-region based spectral analysis is supported as 

well. 
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